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Abstract 

We discuss a kind of generalized concurrence for a class of high dimen- 
sional quantum pure states such that the entanglement of formation is 
a monotonically increasing convex function of the generalized concur- 
rence. An analytical expression of the entanglement of formation for a 
class of high dimensional quantum mixed states is obtained. 

PACS numbers: 03.65.Bz, 89.70. +c 

Key words: Entanglement of formation, Generalized concurrence 

Characterized by the simultaneous occurrence of superposition and correlation in com- 
posite systems, quantum entanglement has in recent years emerged as the key resource in 
quantum information processing [1] and resulted in the explosion of interest in quantum 
computing and communication [2]. To quantify entanglement, a number of entanglement 
measures has been proposed for bipartite systems. Among them, the pioneering contri- 
butions of Bennett et al. [3, 4] defined entanglements of formation and distillation on 
considerations of convertibility vis-a-vis maximally entangled pairs. The relative entropy of 
entanglement [5] is based on distinguishability from the set of separable states. Recently 
there are also discussions of entanglement measure on multipartite case [6]. 

The entanglement of formation is defined for arbitrary dimensional bipartite systems. 

Due to the extremizations involved in the calculation, so far no explicit analytic formulae 
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for entanglement of formation have been found for systems larger than a pair of qubits, 
except for some special symmetric states [7]. 

Let H be an iV-dimensional complex Hilbert space with orthonormal basis e^, % — 1, N. 
A pure state on TC <g> Ti is generally of the form, 

N 

1^) = a iJ e i ® e i> a *i G c (!) 

with normalization X^j=i a ij a ij — 1- The entanglement of formation E is defined to be the 
entropy of either of the two sub-Hilbert space H^H [4] , 

= -Tr ( Pl log 2 Pl ) = -Tr (p 2 log 2 p 2 ) , (2) 

where p\ (resp. p^) is the partial trace of over the second (resp. first) Hilbert space 

ofH^H. Let A denote the matrix with entries given by a^- in (1). p\ can be expressed as 
pi = AA^ . A general density matrix p on H <g> H has pure-state decompositions of \ip a ) of 
the form (1) with probabilities p a , 

M M 

P = J2Pa\^a)(^a\, ^Pa=l (3) 
a=l a=l 

for some M e IN. The entanglement of formation for the mixed state p is defined as the 
average entanglement of the pure states of the decomposition, minimized over all possible 
decompositions of p, 

M 

E(p) = mm J2PaE(\ij a )). (4) 

a=l 

For the case N = 2, (2) can be written as 

where h(x) = — x\og 2 x — (1 — x) log 2 (l — x), C is called concurrence, C{\ip)) = 2\a\\a,2i — 
Oi2«2i|- It is easily verified that E is a monotonically increasing function of C, and hence 
C can be also taken as a kind of measure of the entanglement. Calculating (4) is reduced 
to calculate the corresponding minimum of C(p) = mm J2^f=iPaC(\tp a )), and an analytical 
expression of (4) is obtained [8]. 

For iV > 3, there is no explicit analytical expression for the entanglement of formation in 
general. The concurrences discussed in [9, 10] can be only used to judge whether a pure state 
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is separable (or maximally entangled) or not [11, 12], since the entanglement of formation 
is no longer a monotonically increasing function of these concurrences. 

Nevertheless, for a special class of quantum states, certain quantities (generalized concur- 
rence) were found to simplify the calculation of the corresponding entanglement of formation 
[13]. Namely, if AA^ has only two non-zero eigenvalues Ai and A 2 , each with degeneracy m, 
then 

E(\tp)) = m (-x\og 2 x - (— - x) log 2 ( — - x) \ , (5) 

where 



X = \ \- + \—2^- d2 ) 

2 Im V m 

and d = 2mv / AiA 2 is a kind of generalized concurrence taking values from to 1. One easily 
shows that E(\ip)) is a monotonically increasing function of d. Moreover, E(d) is a convex 
function, 

Q2 E log - 2^T^d? 

w = f|^4 >"■ 

From the monotonicity and convexity the entanglement of formation for a class of high 
dimensional mixed states has been calculated analytically [13, 14] 

In the following we generalize the results in [13] to the case that AA^ has n > 3 different 
non-zero eigenvalues. We present the conditions allowing to derive an explicit lower bound 
of the entanglement of formation for such kind of arbitrary dimensional mixed states and 
calculate the lower bound. 

Let Ai, A 2 , • • • , A„, each with degeneracy m, van < N, be the non-zero eigenvalues of 
AA^ . Aj = \i(u, v), i — 1, 2, . . . , n, are differentiable functions of two real variables u and v. 
We define D = mn\/\i\ 2 • • • A n to be the generalized concurrence. 

[Lemma]. If Aj = \i(u, v), i — 1, 2, . . . , n, satisfy the following conditions: 

E < 0. (6) 

then D is a measure of entanglement in the sense that the entanglement of formation of the 
corresponding pure state is a monotonically increasing function of D. 

[Proof] The normalization condition of Tr(AA^) = = 1, Aj G (0,—), gives 
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rise to 



£-§ = 0. (7) 



The entanglement of formation of is given by 

^(l^)) = -E^ 1 °g 2 A i . (8) 

i 

From (7) and (8) we have 

W = -^ ml092Xi W' (9) 

which is positive if the condition (6) is satisfied. Therefore E{\tp)) is a monotonically 
increasing function of D. ■ 

The Lemma defines a class of pure states, for which a generalized concurrence can be 
still well defined. Here we have supposed that all \ are functions of two real variables. 
This implies that there is only one independent variable, accounting to the normalization 
condition, and the condition (6) could be plausibly satisfied. The most simple case is that 
A x = u, A 2 = v, each with degeneracy to, the generalized concurrence is given by D = 
2mv / wt>, which is just the case discussed in [13], where E{\ip)) is not only a monotonically 
increasing but also a convex function of D, > 0. 

As another example we consider non-zero eigenvalues of AA^ such that Ai = u, A 2 = u+v, 
A 3 = u + 2v, each with degeneracy to, u and v G M taking values (0, ^). The generalized 



concurrence is given by D = 3myu(u + v)(u + 2v). It is straightforward to verify that E is 
a monotonically increasing function of D, since 

v^^Aj, , 1 „ , 1 — 3mv 

E 7^ l °9^ = ~ 75= 1 - 9mV < 0. 

i oD 3to^V3to 1 + 3mv 

Due to the relation 

?£ ( !zP + ^ nAj = 2^ (6m " +/n TrlTO^ )<0 ' 

i? is also a convex function of D. 

As E'd^)) is a monotonically increasing and convex function of D, instead of calcu- 
lating E(p), one may calculate the minimum decomposition (in the sense of (4)), D(p) = 
mm J2^!=i PaD(\4>a)), to simplify the calculations, as long as p has all decompositions on pure 
states with their eigenvalues of AA^ satisfying (6) in Lemma. Nevertheless, like E(\tp)), gen- 
erally the expression of D(\ip)) = to77,VAiA 2 • • • A n could be still quite complicated. 



In fact D{\ip)) is an invariant under local unitary transformations. Associated with a 
general pure state \vfj) given in (1), the following quantities are invariants under local unitary 
transformations [10, 11]: 

N N 

a..-~.n* n ■ _/,, 

JP ' 



Jo = Tr(AA^) = £ dij^j , h = Tr[{AA^f\ = £ a ip a* q a jq a* 



The generalized concurrence defined in [10], 



N " 2 



where C2 = C, does have a very simple form. However, the entanglement of formation is 
generally not a monotonically increasing function of Cat, and Cn can be only used to judge 
wether a state is separable or maximally entangled for the case N > 3. If the generalized 
concurrence D = mny/ X1X2 • • • A n satisfies the Lemma and can be further expressed as 



D — ^fy-fo — h, the calculation of the corresponding entanglement of formation would be 
greatly simplified. 

Let * denote the set of all pure states of the form (1) such that i) the Lemma is satisfied; 
ii) the entanglement of formation is a convex function of D, i.e., X) Y^Wd^ ^W^^ 12 ^ < ^' 



D = mnVAiAa • • • A„ = t%JIo ~ h, that is, T7 /»(«, u) = ^[(m-l)+m 2 £ «)]■ 
We call a mixed state p given by (3) .D-computable if all the decompositions of p on pure 
states belonging to ^ . 

Due to the conditions i) and ii), for a D-computable state p, calculating E(p) is then 
reduced to the calculation of the corresponding minimum of D(p) = min J2^ = iP a D(\i^ a )), 
which simplifies the calculation if D(\ip a )) has a simpler expression than E(\ip a )). The 
condition iii) guarantees that D is a quadratic form of the entries of the matrix A and can 
be expressed in the form of D = \(if)\S'ip*)\ in terms of a suitable matrix S, which allows us 
to find an explicit analytical expression of the entanglement of formation in a way similar 
to the one used in [8] and [13]. It generalizes the results in [13, 14] where the case that AA^ 
has two non-zero eigenvalues is considered. 

Let S ip i q be a symmetric iV 2 x N 2 matrix whose elements are all zero except for 

Sp+N(i-l),q+N(j-l) = S q+N (j-i) jP+N (i-i) = 1, 
Sq+N(i-l),p+N(j-l) = Sp+ N (j-i) jq+N (i-i) = — 1, 



where i,j,p,q = 1,...,N. Let A ipM , A 1 ^ 39 , A % £ 39 and Af 39 , in decreasing order, be the 



eigenvalues of the rank four Hermitian matrix ^ y/pS ip i q p* S ip i q y/p. 

[Theorem]. For a D-computable state p, the minimum decomposition of the general- 
ized concurrence D(p), i.e. the average generalized concurrence of the pure states of the 
decomposition, minimized over all decompositions of p, is given by 

N 

™[ (KY jq - AT q - Af q - Af q f]^ ■ (11) 

«J>P><3=1 

[Proof]. Let r be the rank of p and \vk), k — 1, ...,r, be a complete set of orthogonal 
eigenvectors corresponding to the nonzero eigenvalues of p, such that (vk\vk) is equal to 
the kth. eigenvalue. Other decomposition of p can then be obtained through unitary 

transfer mat ions : 

K> = (12) 
i=i 

where U is a t x t unitary matrix, t > r. We have (wk\S ipjq wf) = (UT ipjq U T )ku where the 
matrix r wjq is defined by T l k P3q = (v k \S ipjq vf). As the matrix S ipjq is symmetric, r tpjq is also 
symmetric and can always be diagonalized by a unitary matrix U such that Ut IP39 U t = 
diag(A l i' )q , ...,A tpjq ) [15]. The diagonal elements A tpjq , a = 1, r, in decreasing order, can 
always be made to be real and non-negative. They are also the eigenvalues of the Hermitian 
matrix R = ^ y /pS^ q p*S i ^ q ^fp. The matrix S ipjq has N 2 — 4 rows and N 2 - 4 columns 
that are identically zero. It can be seen that the corresponding rows and columns of matrix 
S ip i q p*S ipjq are identically zero as well. Thus the N 2 x N 2 Hermitian matrix S ipjq p*S ipjq 



has a rank no greater than four. It follows that the Hermitian matrix y yJJ>S ip ^ q p* S vp3q ^/p 
is ranked at most four: at least Af 39 , A l Q jq , ...,A tpjq are zero. 

There always exits a decomposition consisting of states \wk), k = 1, . . . , r, such that 
(wk\S ipjq Wi) = A^^Ski- Set \yi) = \wi), \yi) = i\w{) for I = 2, ...,r. Any decomposition can 
be written in terms of the states \yk) via the equation \zk) = Z)[=i where V is a t x r 

matrix whose r columns are orthonormal vectors. 

Denote d ip j q = \(ip\S ipjq il)*)\ = 2\a ip aj q — a iq aj P \. The average d ip j q of a general decom- 
position is given by 

(d iPjg ) = E \(vYv T ) kk \ = E|E(^«) 2y «|. ( 13 ) 

k k I 
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where Y is the real diagonal matrix defined by Y ki = (yk\S ipjq yl) ■ Using the fact that 
Efe |(Vfe/) 2 | = 1, one gets 



(d ipjq ) > | T,(Vm) 2 Yu\ > A^ q - A 2 ipjq - As ipjq - A 4 



ipj'i 



ki 



Since 



_ . . , . . mn 



N 



mn 



El 12 
\ a ip a jq ~ a iq a jp\ = 



N 



V d 2 ■ 



using the Cauchy-Schwarz inequality, 

M N N M 

E( E \ms™ q m\ 2 )*>\ E (El(^l^lC)l) 2 ]i 



a=l i,j,p,q=l 



we get 



iV 



Therefore the minimum decomposition of the generalized concurrence D(p) is given by 



A? 



""" [ (A^ 9 — A 2 P3q — A 1 ™ 9 — A t [ 3q ) 2 ]^ 
^ i,j,p,q=i 

Due to convex relation between E{\^)) and D(\ip)), the entanglement of formation of p is 
given by E(D(p)). ■ 

As a simple example we consider a class of pure states on 3 x 3 (which is not the case 
in [13, 14]) with the matrix A given by 

/ an a \2 a i3 \ 



A = 



0,21 &22 a 23 
«21 «22 ^23 



(14) 



/ 



The matrix AA* has two non-zero eigenvalues Ai and A2 satisfying 

AiA 2 = 2(\a u a 2 2 - ai2fl2i| 2 + 1 011*123 - ai3fl2i| 2 + |ai2«23 - «i3fl22| 2 )- 

It is directly verified that all pure states given by (14) belong to ^. 

For all mixed states p with decompositions on pure states (14) (it is easily seen that 
once p has a decomposition with all the pure states given by (14), then all other kinds 
of decompositions belong to ^ too), from the theorem the lower bound of the generalized 
concurrence D(p) is given by 

v^KA^-A^-A^-Al 122 )^^ 1 ^ 



We have studied the entanglement of formation for higher dimensional quantum mixed 
states. It has been shown that under certain conditions the entanglement of formation 
is a monotonically increasing and convex function of a generalized concurrence. For a 
class of (D-computable) arbitrary dimensional mixed states an explicit lower bound of the 
entanglement of formation is derived. The generalized concurrence defined in this note is 
a generalization of the one in [13]. However, as we imposed on the condition Y\fi(u,v) = 

i 

^[(m — 1) + m 2 fi(u,v)fj(u,v)] in calculating the lower bound of the entanglement of 

i¥=3 

formation, the .D-computable states here only cover a part of the d-computable states in 
[14]. In fact, in order to get the lower bound of the concurrence by using the method above, 
it is not necessary for D to be of the form Y^j p q =\ \ a ip a jq ~ a-i q Q-jp\ 2 ^ but a form of \F (a)| 2 
would be enough, where F(a) stands for an arbitrary quadratic polynomial of the elements 
dij . 

We have assumed that Aj ^ 0, % — 1, ...n, in our theorem. For n > 1, all the pure states 
in \1> are entangled. Therefore the D-computable states are not separable in this 
all possible decompositions, including the minimum decomposition of the entanglement of 
formation are in ^. One can check if these states are bound entangled by simply checking 
whether their are PPT (positive partial transposition) or not. 

References 

[1] C. H. Bennet, Phys. Scr. T76, 210 (1998). 

[2] M. Nielsen and I. Chuang, Quantum Computation and Quantum Information (Cam- 
bridge University Press, Cambridge, England, 2000). 

[3] C. H. Bennett, G. Brassard, S. Popescu, B. Schumacher, J. A. Smolin, and W. 
K.Wootters, Phys. Rev. Lett. 76, 722 (1996); 

C. H. Bennett, D. P. DiVincenzo, J. A. Smolin, and W. K. Wootters, Phys. Rev. A 54, 
3824 (1996). 

[4] C.H. Bennett, H.J. Bernstein, S. Popescu, and B. Schumacher, Phys. Rev. A 53, 2046 
(1996). 

8 



[5] V.Vedral and M. B. Plenio, Phys. Rev. A 57, 1619 (1998), and references therein. 
V. Vedral, Rev. Mod. Phys. 74, 197 ( 2002). 

[6] M.H. Partovi, Phys. Rev. Lett. 92, 077904 (2004). 

[7] B.M. Terhal, K. Gerd and K.G.H. Vollbrecht, Phys. Rev. Lett. 85, 2625 (2000). 

[8] S. Hill and W.K. Wootters, Phys. Rev. Lett. 78, 5022 (1997). 
W.K. Wootters, Phys. Rev. Lett. 80, 2245 (1998). 

[9] A.Uhlmann, Phys. Rev. A 62, 032307 (2000). 

P. Rungta, V. Buzek, CM. Caves, M. Hillery, G.J. Milburn, Phys. Rev. A 64, (042315) 
(2001). 

[10] S. Albererio and S.M. Fei, J. Opt. B: Quantum Semiclass. Opt. 3, 1-5(2001). 

[11] S. Albeverio, S.M. Fei and D. Goswami, Phys. Lett. A, 91-96 (2001). 

[12] S.M. Fei, X.H. Gao, X.H. Wang, Z.X. Wang and K. Wu, Phys. Lett. A 300, 559- 
566(2002). 

[13] S.M. Fei, J. Jost, X.Q. Li-Jost and G.F. Wang, Phys. Lett. A 310, 333-338(2003). 
[14] S.M. Fei and X.Q. Li-Jost, Rep. Math. Phys. 53, 195-210(2004). 

[15] R.A. Horn and C.R. Johnson, Matrix Analysis, Cambridge University Press, New York, 
1985. 



9 



